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ABSTRACT 

This paper presents the results of clock synchronization experiments 
using OMEGA transmissions from North Dakota on IS.lOkHz and 
12 . 85 kHz . The OMEGA transmissions were monitored during April 
1974 from NASA tracking sites located at Madrid, Spain; Canary Island; 
and Winkfield, Finland. The sites are located at distances between 
6600 kilometers (22, 100 us) to 7300 kilometers (24,400^8) from North 
Dakota. 

The data shows that cycle identification of the received signals was ac- 
complished. There are, however, discrepancies between the measured 
and calculated propagation delay values which have not been explained, 
but seem to increase with distance between the receiver and the trans- 
mitter. The data also indicates that three strategically located OMEGA 
transmitting stations may be adequate to provide worldwide coverage 
for clock synchronization to within ± two (2) microseconds . 


INTRODUCTION 

The field tests were conducted from April 1 through April 22, 1974. The ob- 
jectives of these tests were to determine the adequacy of unique OMEGA trans- 
missions for time determination at distances greater than 6000 km from the 
transmitter and to collect experimental data that would indicate the minimum 
number of sites, appropriately located, required for worldwide clock synchro- 
nization to within ±2 us . 

Receiving transmissions from the OMEGA North Dakota (N. D.) station, data were 
collected at Madrid, Spain (7185 km); Canary Island (7298 km); and Winkfield, Eng- 
land (6626 km). Previous test results conducted at Greenbelt, Md. (1934km); 
Washington, D.C. (1922km) and Rosman, N.C. (1794km) were reported at the 1973 
PTTI Meeting Further tests will be conducted during September-October 1975 
at Hawaii and Australia. Transmissions from OMEGA stations N.D. andHawaii 
will be used. Transmissions from OMEGA N.D, will be monitored at Kauai, Ha- 
waii (5926km) and from both N.D. and Hawaii will be monitored at Orroral, Aus- 
tralia (14,408km and 8443km respectively from the transmitters). 
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The Use of the OMEGA system for time transmission not only will augment ex- 
isting systems such as Loran-C by providing an additional capability in the north- 
ern hemisphere, but would also provide coverage in the southern hemisphere 
where Loran-C transmitters do not exist and time synchronization is limited to 
within ±25 /zs. 

PROPAGATION DELAY MEASUREMENT 


If the phase velocities for two VLF signals (f^ and fj ) are v^ and Vj respectively, 
then the propagation delays for a path length D are: 


and 




( 1 ) 

( 2 ) 


Since the two signals were emitted in phase at time t^^^j at the transmitter , the 
signals as received at a station at time t^^ relative to the station clock must sat- 
isfy the following conditions'^): 


txm - trvi = tp™ = n + (nj + Anj ) Ti (3) 

txm~frv2 ~ “ UTg ±(n2 +An2)T2 (4) 

where t^^ and tjy^ are the time of s^al reception of f j and f^ relative to the 
receiving station clock, and tp® and tp™ are measured propagation delay values. 
The received cycle of each carrier, n^ and n 2 is determined from the phase 
measurements Anj and An 2 - 

The ambigui^ identification of n for (20 milliseconds) is considered a priori 
knowledge or determined from the theoretical predicted propagation delay (see 
table 1) which is in error by much less than 4 milliseconds. The receiver de- 
sign actually requires the ambiguity resolution to 4 milliseconds (see cycle iden- 
tification p. 193 and reference 2). 


Equations (3) and (4) can be used to calculate the measured propagation delay if 
the receiving station clock time is known relative to the transmitter clock. Hav- 
ing measured the propagation delays via portable clock, the received VLF sig- 
nals can be used to synchronize a receiving station clock. 
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Table 1 


Ambiguity Identification of n from Predicted Propagation Delays 


Ambiguity Determination 


NASA 

Tracking 

Sites 


Madrid 


Winkfield 


Predicted Delay (ms) 

tpi 

(13.10 kHz) 

*P2 

(12.85 kHz) 

23,974.8 

1 

23,972.1 

24,346.3 

24,343.5 

22,113.0 

22,110.5 






nTg (MS) 


20,000 


20,000 


20,000 
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CYCUES CYCLES 



MEASURED AND PREDICTED DIURNAL PHME OP 
13. 1 KHz OMEGA H DAKOTA TO MADRID , SPAIN 



MEASURED AND PREDICTED DIURNAL PHASE OF 
I2.85KHZ OMEGA N- DAKOTA TO MADRID. SPAIN 


Figure 2 . Diurnal phase records of OMEGA 
North Dakota transmitted signals as received at 
Madrid, Spain on April 3 to 7, 1974. 



MEASURED AND PREDICTED DIURNAL PHASE OF 
I3.I KH1 OMEGA N. DAKOTA TO CANARY ISLANO, SPAIN 



MEASURED AND PREDICTED DIURNAL PHASE OF 
12.85 KHz OMEGA N. DAKOTA 10 CANARY ISLANO. SPAIN 

Figure 3. Diurnal phase records of OMEGA North 
Dakota transmitted signals as received at Canaiy 
Island , Spain on April 6 to 11, 1974. 
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MEASURED AND PREDICTED DIURNAL PHASE OF 
IS. I KHz OMEGA N. DAKOTA TO WINKFIELD. ENGLAND 



MEASURED AND PREOICTEO DIURNAL PHASE OF 
12.85 KHz OMEGA N. DAKOTA TO WINKFIELD . EI^LAND 

Figure 4. Diurnal phase records of OMEGA North 
Dakota transmitted signals as received at Winkfield, 

England on April 16 to 20, 1974. 

The digital data were grouped in time periods according to the diurnal phase rec- 
ord. This was done by taking the average of tiie data points obtained over 15 or 
20 minute intervals during the periods of sunrise, daytime, sunset, and nighttime 
(see table 2). Due to local interference problem in Madrid, only phase record 
during the da 3 rtime transmission period were collected. 

With reference to the results given in table 3, An^ is the phase difference of the 
received signal (13.10 kHz) relative to the station clock and similarly for Au 2 
(12.85kHz). The difference of the two received signal phases is given as An^j 
which is AOj - . 


CYCLE DETERMINATION 
Combining equations (3) and (4) : 

Atp s At® - t® = (ni + Ani)Ti -(n2 +An2)r2' (5) 


Let: 


m = Uj - nj 
At - Ti- t2 
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Table 2 


Time of Day (GMT) When Data Were Collected 


Time of Day 

Madrid* 

Canary Island** 

Winkfield** 

Sunrise 

— 

0630 - 1050 

0510 - 1150 

Day 

1100-1800 

1110 - 1810 

1210 - 1650 

Sunset 

— 

1830 - 0250 

1710 - 2350 

Night 


0310 - 0610 

0010 - 0450 


*Data was collected at IS minute intervals. 
**Data was collected at 20 minute intervals. 


Table 3 

Average of Measured VLF Signal Phase in April 1974 


Averting 

Period 

Madrid 
(April 4-7) 

Canary Island 
(April 6-11) 

Winkfield 
(April 16-19) 

Anj 


An^ 

AUi2 

AUj 

Anj2 

Sunrise 

— 

— 

— 

.140 

— 

.557 

Day 

.825 

.899 

.063 

.141 

.508 

.556 

Sunset 



— 

.136 


.554 

Night 

— 

i 

.808 

.164 

.999 

.556 


and 


_ - n T2 _ Ti T2 

® T2 - Tj At 


Anj2 = AUj - An^ 


then equation (5) becomes: 
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(6) 


At“ = (m + Anj ) (t2 + At) - (n2 + An2 ) T 2 

= (n^ - 112 + Ahj - A112 ) T 2 + (ill + Ani )At 

At^ = (m + Anj2)T2 + (ni +Aiii)At 


Rewriting equation (6); 


T 2 

ni +Ani = - (m + Ani2) ^ 

Th 

= (m + Ani2) ^ - rT=^ 

Similarly: 

Th At? 

n2 +An2 = (mtAni2) ^ (8) 

Either equation (7) or (8) can be used to determine the received carrier cycle of 
either or £2 respectively. This is so because Aui, An 2 and Ani 2 are measured 
quantities; m is the highest integer miiltiple of the beat frequency period within 
the propagation delay value between the transmitter and the receiving station. 
Values of (n) and (m) are therefore known. The values of (n) and (m) used for 
the stations visited in 1973 and 1974 are given in table 4. 

The particular cycle of a carrier signal received at a station as calculated from equa- 
tion (7) or (8) is found by assuming a parametric value for At™, i. e. , At™ = 0, 

1, 2, 3, 4, etc. iiB. After n^ + Ani has been calculated, the value is rounded to 
an integer by the following rule: 

(a) Round n^ + An^ downward to n^, if 

Auj < 0.44546 -»• 0.4455 0.446 0.45 ^ 0.4 ^ 0 

(b) Roimd ni + An ^ upward to n^ + 1, if 

Anj ^ 0,54546 -»• 0. 5455 0.546 0.55 0. 6 1 

The 0.1 cycle difference, between the upper bound of one cycle, n^ , (0.44546) 
and lower botmd of another cycle, ni + 1, (0,54646), is the wall thickness of the 
cycle well. This is the same as sayitsg that 90% of the data is grouped correctly 
into a particular cycle and 10% of the data is grouped into either one cycle high 
or one cycle low. 
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Table 4 


(n) and (m) Values Used for Cycle Determination 
for OMEGA North Dakota and NASA Tracking Stations 


Station 

Predicted 
Prop^ation 
Delay (ns) 

Ambiguity Reduction 

n 

m 

Predicted 
Proi^Lgation 
Delay <^s) 

Rosman, N.C. 

5994 

0 

1 

1994 

USNO, Wash. , D. C. 

6421 

0 

1 

2421 

GSFC, Md . 

5458 

0 

1 

1458 

NELC, Ca. 

7388 

0 

1 

3388 

Winkfield 

22113 

1 

0 

2113 

Madrid 

23974 

1 

0 

3974 

Canary Island 

24346 

1 

i 

1 

346 


PROPAGATION DELAY ANOMALY, At® 


The knowledge of the propagation delay anomaly in the frequency range of 11. 05 
kHz and 13. 15kHz is not accurate. This is because there is a lack of adequate 
phase velocity data in this frequency range. A more accurate method to deter- 
mine the propagation delay anomaly is perhaps by the use of the time measure- 
ment. Table 5 gives a comparison of predicted and measured propagation delay 
values . 

One area in which more research needs to be done is the determination or anal- 
ysis of propagation anomalies in the OMEGA frequency range where propagation 
predictions are not known to sufficient accuracy for microsecond clock synchro- 
nization. The propagation delays must be measured with a clock which is known 
relative to the clock at the transmitter . An alternate approach is to obtain an 
average of the measured phase differences for fi and f 2 for day and n^ht time 
transmissions . From these average values one can calculate the propagation 
delay, such as shown by Table 6. Table 6 gives the sample calculation of cycle 
determination for data obtained from Winkfield. 
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Table 5 


Comparison of Predicted and Measured Propagation Delay Values in 
Microseconds for Indicated Values of tp’s 


Station 

Predicted 

Measured 

Received Carrier (f j ) 
C^cle Determination 

*Pi 

Atp 

n 

m 

Mp 


Predicted 

Measured 

Winkfield 

22,113.0 

2.5 

1 

0 

0 

22,253.7 

— 

— 






3 

22,099.5 

27 

27 

Madrid 

23,974.8 

2.7 

1 

0 

0 

23,644.8 








3 

23,498.1 

52 

45 

Canary 

24,346.3 

2.8 

1 

1 

0 

24,585.7 



Island 





3 

24,508.6 

57 

59 






4 

24,432.3 

— 

58 

i 





5 

24,355.9 

— 

57 

Greenbelt 
i 

6,458 

— ! 

0 

1 

2 

6,463.8 

84 

84 


RESULTS 

The propagation delays for each station visited were measured. The results are 
given in Table 5. It is to be noted that for Atp «= 3 ms, the measured propagation 
delays for Winkfield and Madrid are lower than the predicted values; Canary 
Island is higher. 

A sample calculation of cycle determination of n^ for Atp == 0, 1, 2, 3, 4 ms is 
given in Table 6. The average of the measured propagation delays are 2097.5 
and 2137.3 /us respectively fof day'^ and night-time transmission paths. 

It is noted that the propagation anomaly for a particular receiving station is de- 
termined only when the station clock is known relative to the transmitter clock. 
This information can be obtained by the use of a portable clock. Ourii^ these 
field tests the OMEGA North Dakota station was measured relative to the U.S. 
Naval Observatory (USNO) Master Clock (MC) via a portable clock. It was de- 
termined by the USNO that OMEGA Noith Dakota was slow relative to USNO-MC 
by 7.4 microseconds. Also, during the tests a portable clock referenced to the 
USNO-MC was used. All known biases such as clock differences must be taken 
into accoimt in the calculation of the propagation delay values. Once these values 
are established, microsecond time can be obtained since the only limitation is 
the phase stability of the received VLF signals . 
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Tabl e 6 


Sample Calculation of Cycle Determination for data 
Obtained from Winkfield, England » 


Date 

Data 

An^ 


nj For Atp = 

(MS) = 

No 

Cycle 

Cycle 

0 

1 

2 

3 

4 

(27 + Anj) Tj 

1974 Apr 

16 

9 

0.487 

0.562 

29 

29 

28 

27 

27 

2099.2 

1310 - 1750 

17 

13 

0. 510 

0.548 

29 

28 

27 

27 

26 

2100.0 

(Day Time 

18 

15 

0.508 

0.556 

29 

28 

28 

27 

26 

2099.8 

Path) 

19 

15 

0. 498 

0. 560 

29 

29 

28 

27 

27 

2099.1 

1974 Apr 

16 

15 

1.088 

0. 556 

29 

28 

28 

27 

26 1 

2144.1 

0110 - 0530 

17 

15 

1.003 

0.554 

29 

28 

28 

27 

26 

2137.6 

(Night Time 

18 

15 

0. 955 

0.566 

30 

29 

28 

28 

27 

2133. 9 

Path) 

19 

15 

0.950 

0.547 

29 

28 

27 

27 

26 

1 

2133.6 


t“ (Day) - (27 + Anj ) >1 = 2097. 5 ms 


t® (Night) - (27 + Anj) = 2137. 3 ms 

CONCLUSIONS 

]^esults based on the phase data collected during the field tests show that cycle 
determination of a received carrier signal can be determined and that the accu- 
racy of time measurement using two VLF signals is within ±2 microseconds . 

For a worldwide time trans m ission system the field test results Indicate that a 
•minimum of Ihree transmission stations are needed. These stations which trans- 
mit the present dual frequency format must be selected from the OMEGA Navi- 
gation System appropriately in order to provide the coverage. 
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